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ABSTRACT 


The  H-field  in  textbooks  is  derived  by  a  finite  current  element  along  the  z-axis.  This 
technical  memorandum  derives  the  H-field  from  a  current  element  at  arbitrary  orientation  in 
terms  of  the  endpoints  of  the  current  element.  This  result  is  useful  for  a  revision  of  IEMCADS  as 
part  of  the  EMC  effort. 
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MAGNETIC  FIELD  INTENSITY  DUE  TO  A  LINE  OF  CURRENT  OF  ARBITRARY 
ORIENTATION  IN  CARTESIAN  COORDINATES 


Textbooks  give  the  magnetic  field  intensity  of  a  line  of  current  in  cylindrical  coordinates  as 


H  = 


4^  (sina2  -  sinai)  a<j, 


(1) 


where  the  line  of  current  is  along  the  z-axis .  The  angles  CX2  and  ai  are  given  in  figure  1.  The 
quantity  p  is  perpendicular  distance  to  the  z-axis.  The  angular  unit  vector  in  cylindrical  coordinates 
is  a<j).  However,  for  use  in  the  Intelligent  EMC  Analysis  and  Design  System  (IEMCADS),  the 
magnetic  field  intensity  has  more  utility  if  expressed  in  cartesian  coordinates  and  alloting  for 
arbitrary  orientations  of  the  line  of  current. 

Figure  2  gives  the  configuration  of  the  line  current  and  coordinates  of  the  final  result.  The 
final  result  gives  the  magnetic  field  intensity  in  terms  of  the  endpoints  of  the  line  current,  the 
current,  and  the  perpendicular  distance  to  the  axis  of  the  line  current.  The  derivation  of  this  result 
is  as  follows: 

The  vector  ft,  has  a  direction  given  by 

ft,  =  (X2  -  Xi)  ax  +  (Y2  -Yi)ay  +  (Z2  -  Zj)az  .  (2) 

The  final  expression  must  be  in  terms  of  the  distance  d  perpendicular  to  the  line  of  charge.  To  do 
this,  the  point  is  introduced  so  that  a  line  drawn  from  the  point  where  the  field  is  completed  or 
measured  is  perpendicular  to  the  axis  of  the  line  of  current.  In  figure  2  this  means  we  introduce  the 
point  P3(X3,Y3,Z3)  along  the  axis  of  the  line  current  such  that  the  vector  P3P1  is  perpendicular  to 
P3P4  =  d  where  4  is  the  point  of  field  computation  or  measurement.  These  vectors  are  given  by 

P3P4  =  (X4  -  X3)  aX  +  (Y4  -Y3)ay  +  (Z4  -  Z3)az  (3) 

P3P1  =  (Xi  -  X3)  aX  +  (Yl  -Y3)aY  +  (Zi  -  Z3)az .  (4) 

The  mathematical  requirements  for  the  location  of  point  3  are  given  by  the  dot  product 

P3P1  •  P3P4  =  0  (5) 

and  the  equation  of  the  straight  line  defining  this  axis  is 

X3-X1  Y3-Y1Z3-Z1  (6) 

a  b  c 

Equation  (6)  c.an  be  written  as  two  equations 

X3  =  l  (Z3-ZO  +  X1  (7) 


and 


1 


(8) 


Y3  =  ~  (Z3-ZD  +  Y1. 


Equation  (5)  expressed  explicitly  in  cartesian  coordinates  is 

(X4  -  X3XX1  -  X3)  +  (Y4  -  Y3XY1 .  Y3)  +  (Z4  -  Z3)(Zi  -  Z3  )  =  0 
Substituting  7  and  8  into  9  yields 

[x4-|f  (Z3  -Zi)  +  X ij  ]  [- | (Z3 -Zi)] 

+  [x4  -  {£  (Z3  -  Zi)  +Y 1}  ]  [-  £  (Z3  -Zi)] 


(9) 


+  [Z4-Z3][Zi-Z3]=0 


(10) 


This  equation  is  quadratic  in  the  unknown  Z3  and  can  be  solved  by  the  quadratic  equation  since  X\, 
Yi,  Z,  Xi.  Yi,  Zi,  andXq,  Y4,  Z4  are  considered  known  quantities.  Once  equation  (10)  is  solved 
for  Z3,  X3,  and  Y3  are  easily  obtained  from  equations  (7)  and  (8). 

The  solution  of  X3,  Y3  and  Z3  combined  with  the  known  solution  of  X4,  Y4,  and  Z4  will 
yield  the  perpendicular  distance  d  from  the  point  of  field  measurement  to  the  axis  of  the  line 
current.  The  quantity  d  is  given  by 


d  =  V(X3  -  X4)2  +  (Y3  -  Y4)2  +  (Z3  -  Z4)2  .  (11) 

The  magnetic  field  intensity  for  the  arbitrary  orientation  in  cartesian  coordinates  is  finally 

H  =  [sinp2a±2  -  sinpian]  ( 12) 

4rcd 


where 

as) 

and 

sinPi=afi-  (14) 

The  distances  d32,  d24,  d3i,  and  di4  are  given  by 

d32  =  V(X3  -  X2)2  +  (Y3  -  Y2)2  +  (Z3  -  Z2)2  (15) 

d31  =  V(X3  -  Xi)2  +  (Y3  -  Yi)2  +  (Z3  -  Zi)2  (16) 

d24  =  V(Xl  -  X4)2  +  (Yi  -  Y4)2  +  (Zi  -  Z4)2  (17) 
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dl4  =  V(Xl  -  X4)2  +  (Y 1  -  Y4)2  +  (Zi  -  Z4)2 
The  unit  vectors  a_L  are  given  by 

p3P2  x  P3P4 

a±2-  |p3p2  x  p3p4| 


(18) 


(19) 


and 


a,  .  p3Pl x  P3P4 
a±1  “  IP3P1  x  P3P4I  ' 


(20) 


These  unit  vectors  can  explicitly  be  written  in  terms  of  the  end  points  of  the  current 
element  and  the  now  known  coordinate  X3,  Y3,  and  Z3. 


To  do  this  we  express  P3P1,  P3P4,  and  P3P2in  known  cartesian  coordinates 
P3P1  =  (X3  -  Xi)&x  +  (Y3  -  Yi)ay  +  (Z3  -  Zi)fiz  (21) 

P3P4  =  (X3  -  X4)aX  +  (Y3  -  Y4)ay  +  (Z3  -  Z4)az  ‘  (22) 

P3P2  =  (X3  -  X2)ax  +  (Y3  -  Y2)ay  +  (Z3  -  Z2)aZ  .  (23) 


Taking  the  cross  products  one  merely  takes  the  determinant  of  the  following  matrices 

P3P1  x  P3P4  = 


ax 

(X3-Xi) 
(X3  -  X4) 


ay 

(Y3-Yj) 

(Y3-Y4) 


P3P2  X  P3P4  = 


ax 

(X3  -  X2) 
(X3-X4) 


ay 

(Y3-Y2) 

(y3-y4) 


aZ 

<z3  -  Zi) 
(Z3  -  Z4) 


Taking  the  determinant  of  equation  (24)  yields  (25) 
=  aX[(Y3  -  Yi)(Z3  -  Z4)  -  (Y3  -  Y4)(Z3  -  Zi)] 

+  ay[(Z3  -  Zi)(X3  -  X4)  -  (X3  -  XL)(Z3  -  Z4)] 

+  aZ[(X3  -  Xi)(Y3  -  Y4)  -  (Y3  -  Yi)(X3  -  X4)] 


(24) 


aZ 

(Z3  -  Z2) 
(Z3-Z4) 


(25) 


(26) 


Taking  the  determinant  of  equation  (26)  yields 
=  aX[(Y3  -  Y2)(Z3  -  Z4)  -  (Z3  -  Z2)(Y3  -  Y4)] 

+  ay[(Z3  -  Z2)(X3  -  X4)  -  (X3  -  X2)(Z3  -  Z4)] 

+  az[(X3  -  X2)(Y3  -  Y4)  -  (X3  -  X4XZ3  -  Z2)]  (27) 
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Taking  the  absolute  values,  we  obtain 

IP3P1 x  P3P4I 

=  ([(Y3  -  Y1XZ3  -  Z4)  -  (Y3  -  Y4)(Z3  -  Zi)]2 
+  [(Z3  -  Zi)(X3  -  X4)  -  (X3  -  X1XZ3  -  Z4)]2 

+[(X3-XiXY3-Y4)-(Y3-YiXX3-X4]2)1/2  (28) 

IP3P1 X  P3P4I 

=  ([(Y3  -  Y2XZ3  -  Z4)  -  (Z3  -  Z2)(Y3  -Y i)]2 
+  [(Z3  -  Z2XX3  -  X4)  -  (X3  -  X2)(Z3  -  Z4)]2 
+[(X3  -  X2XY3  -  Y4)  -  (X3  -  X4)(Y3  -  Y2]2)l/2 
The  final  form  for  the  H-field  at  arbitrary  orientations  is 

H  =  -7—  [sinp2a±2  -  sinpian] 

4jtd 

where  d  is  given  by  equation  (1 1),  p2  and  Pi  are  given  in  equations  (13)  through  (18),  and 
aj_is  given  by  equations  (19)  through  (29). 

CONCLUSIONS 


(29) 

(30) 


The  final  expression  for  the  H-fields  is  given  by  equation  (30)  and  all  the  appropriate 
substitutions.  The  final  result  gives  the  field  at  some  arbitrary  point  in  terms  of  the  endpoints  of 
the  current  element,  the  current,  and  the  perpendicular  distance  to  the  axis  of  the  current  element. 
This  result  differs  from  the  current  elements  presently  in  IEMCADS  in  that  it  gives  the  resulting 
field  with  the  current  element  at  arbitrary  orientations.  IEMCADS  is  presently  limited  at  giving 
fields  for  current  elements  fixed  along  the  z-axis  and  consequendy  is  of  limited  value.  This  result 
will  be  incorporated  into  IEMCADS  to  give  greater  flexibility  in  its  EMI  predictions. 
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Figure  1.  The  H-field  caused  by  a  finite  current  element  on  the  z-axis. 
The  H-field  at  Point  2  is  given  by  equation  (1). 
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